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Abstract 



The decay K~ — > //"z/tt has been studied using in-flight decays detected with "ISTRA+" setup 
operating in the 25 GeV negative secondary beam of the U-70 PS. About 112K events were 
used for the analysis. The A+ and Aq slope parameters of the decay formfactors f+{t), fo{t) 
have been measured : 
A+ = 0.0321 ± 0.004(stat) ±0.002(syst) 

Ao = 0.0209 ± 0.004(stat) ±0.002(syst); the correlation dXo/dX+ = -0.46 
The limits on the possible tensor and scalar couplings have been derived: 
fr/f+iO) = -0.021 ± 0.028(stat) ±0.014(theory); 
fs/f+{0) = 0.004 ± 0.005(stat) ±0.005(theory) 
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Figure 1: The layout of the "ISTRA+" setup. 



1 Introduction 

The decay K —>■ fiu7r^{K^^) is known to be a key one in hunting for phenomena beyond the 
Standard Model (SM). In particular, significant efforts have been invested into T- violation 
searches by the measurements of the muon transverse polarization err- In our analysis, based 
on ~ 112K events of the decay, we present new search for S and T interactions by fitting the 
Kfj^s Dalitz plot distribution, similar to as it was done for the Ke^ decay: Another subject 
of our study is the measurement of the V-A f+(t) fo(t) formfactor slopes A+ and Aq. 

2 Experimental setup 

The experiment is performed at the IHEP 70 GeV proton synchrotron U-70. The experimental 
setup "ISTRA+" has been described in some details in our recent paper on K^s decay [Q]. 
A schematic view of the detector is shown in Fig.l. The setup is located in the 4A negative 
unseparated secondary beam. The beam momentum is ~ 25 GeV with Ap/p ~ 2%. The 
admixture of K~ in the beam is ~ 3%. The beam intensity is ~ 3 ■ 10^ per 1.9 sec. U-70 spill. 

3 Event selection 

During 3 weeks physics run in March-April 2001, when the muon identification was in full 
operation, 363M events were logged on DLT's. This information is supported by about lOOM 
MC events generated with Geant3 0. Some information on the reconstruction procedure is 
presented in here we touch only points relevant for the i^^^s events selection. 

The muon identification is based on the information from the SPi - a 576-cell lead glass 
calorimeter and HC- a scintillator-iron sampling hadron calorimeter, subdivided into 7 longi- 
tudinal sections 7x7 cells each 0. The calorimeters are located at the very end of the setup, 
after the main magnet (M2) and the last elements of the tracking system: drift tubes (DT) and 
the matrix scintillation hodoscope (MH). The first requirement is that the energy of the SPi 
cluster, associated with the charged track is less than ~ 2.5 MIP's; the HC energy, associated 
with the track should also be less than 2.5 MIP's. The last selection requires that more than 
10% of the HC associated energy is deposited in 2 last layers (out of 7) of the HC. The efficiency 
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Figure 2: The 77 mass spectrum for the events with the identified muon and two extra showers. 

of the algorithm to muons is tested on K ^ ixv events and is found to be ~ 70%. The vr ^ /i 
misindification is measured on K~ — > 7r~7r'' decay and is ~ 3%. After the muon identification, 
the selection of the events with two extra showers results in M^^ spectrum shown in Fig. 2. 
The TT*^ peak has a mass of M^-o = 134.6 MeV, and a resolution of 8.6 MeV. The missing mass 
squared- {Pk — Ppi — P-k^Y-, where P are the corresponding four-momenta, is presented in Fig. 3. 
The cut is ±0.01 GeV^. The further selection is done by the requirement that the event passes 
2C K — > ixvTT^ fit. The missing energy Ek — E^ — E^o after this selection is shown in Fig. 4 The 
peak at low Emiss corresponds to the remaining K~ ti^tt^ background. The corresponding 
cut is Emiss > 1.4 GeV. The surviving background is estimated from MC to be less than 4%. 
The detailed data reduction information is shown in Table. 1. 

4 Analysis 

The event selection described in the previous section results in selected 112K events in 2001 
data. The distribution of the events over the Dalitz plot is shown in Fig. 5. The variables 
y = 2E^/Mk and z = 2Et^/Mx, where E^^, E^^ are the energies of the muon and 7r° in the kaon 
c.m.s are used. The most general Lorentz invariant form of the matrix element for the decay 
K i2U7T° is g: 

M = (1 + 7') [mKfs - I [{Pk + P.)aU + (Pk - P.U-]^ + i^a^pP^P^]v{p,) 

(1) 



X^/ndf 44.89 / 9 



PI -5847. 
P2 0.1345E+06 
P3 -0.5623E+06 
P4 6250. 
P5 0.1346 
J^6 0.8602E-02 
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Figure 3: The missing four-momentum squared {Pk — — P^o)'^ for the selected events . The 
points with errors are the data, the histogram- MC. 
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Figure 4: The missing energy for the fiir^ events. The points with errors are the data, the 
histograms- MC. The dark(blue) peak at zero value corresponds to the MC-predicted K tt^tc^ 
background. The arrow indicates the cut value. 



4 



Table 1: Event reduction statistics. 



Run 


March- April 2001 


Nevents OU tapCS 


363.002.105 


Beam track reconstructed 


268.564.958 =74 % 


One secondary track found 


134.227.095 =37% 


Written to DST 


107.215.783 =30 % 


/i~ identified and 7r° identified 


218.813 


|Ml,J <0.01 


195.799 


K yU,z/7r° accepted 


166.495 


Erniss > 1-4 GeV 


112.157 



It consists of scalar, vector and tensor terms, fs, fr, f± are functions oft= {Pk — P-k)"^- In the 
Standard Model (SM) the W-boson exchange leads to the pure vector term. The "induced" 
scalar and/or tensor terms, due to EW radiative corrections are negligibly small, i.e the nonzero 
scalar/tensor form factors indicate a physics beyond SM. 

The term in the vector part, proportional to /_ is reduced(using the Dirac equation) to a 
scalar formfactor. In the same way, the tensor term is reduced to a mixture of a scalar and a 
vector formfactors. The redefined /+(V), -Fs'(S) and the corresponding Dalitz plot density in 
the kaon rest frame p{Et^,E^) are H]: 



V = f+ + {m^/mK)h 



S = fs + {m,/2mK)f^+{l + ^ ^ - — ) /t 

2mj^ rriK rn^ 



ml_ 2E„ E^ 

.2 



piE^,E^) ~ A-\V\' + B-ReiV*S) + C-\S\' (2) 

A = mK{2E^E,-mKAE^) ~ml{E,-^AE^) 

B = m^mK{2E^ - /\E^) 

C = mlAE^; AE^ = E™"" - K"" = — ^ 



2m 



K 



Following 0] a scalar formfactor /o is introduced: /o(t) = f+{t) + „^i '^_,^2 f-{t) and linear 

dependence of /+, /o on t is assumed: = /+(0)(1 + \+t/ml); /o(t) = /+(0)(1 + Xot/ml). 

Then/_ = /+(0)(Ao-A+)^^^. 

The procedure for the experimental extraction of the parameters A+, Aq, fs, /r starts from 
the subtraction of the MC estimated background from the Dalitz plot of Fig. 4. The background 
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Figure 5: Dalitz plot {y = 2E^/Mk', z = 2Et^o / Mk) for the selected K —>■ fiuir^ events after 
the 2-C fit. 

normalization was determined by the ratio of the real and generated tt^tt^ events. Then 

the Dalitz plot was subdivided into 20 x 20 cells. The background subtracted distribution 
of the numbers of events in the cells (i,j) over the Dalitz plot, for example, in the case of 
simultaneous extraction of A+, Aq and j^^, was fitted with the function: 

p{^,J)- E W,,,,,,{^,J)■\1^■\',^■{fs/Um'' (3) 

ki;ki+k2+k3=0,l,2 

Here W^fcjfc2fc3 are MC-generated functions, which are build up as follows: the MC events are 
generated with constant density over the Dalitz plot and reconstructed with the same program 
as for the real events. Each event carries the weight w determined by the corresponding term 
in formula 2, calculated using the MC-generated("true") values for y and z. The radiative 
corrections according to [0] were taken into account. Then VFfc^fcjfcs is constructed by summing 
up the weights w of the events in the corresponding Dalitz plot cell. This procedure allows to 
avoid the systematic errors due to the "migration" of the events over the Dalitz plot because 
of the finite experimental resolution. 

5 Results 

The results of the fit are summarized in Table. 2. 

The first line corresponds to pure V-A SM fit. The first column is independent fit of our K^^ 
data. The A+ ^ Aq correlation parameter is: ^ = -0.46. The A+ value = 0.0321 ± 0.004 
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Table 2: Results of the fit. 







yU UTT^ + e z/7r° 


A+ 
Ao 


n n'io 1+0-0040 

U-UOZi_o.0040 


0.0296t°:°°ll 
0.0209t°:°°i 


A+ 


n riQ0 1+0-0040 

U.UOZi_o,o040 


0.0297tHS}l 


Ao 


0.01700 


0.01700 


/5//+(0) 


n nn^4+o-oo5e 

U.UUO4_0.0058 


0.0039;°:°°! 


A+ 


n n'?'?s+0-0037 

U-UO-30_o.0037 


o.o299;°:°°ll 


Ao 


0.01700 


0.01700 


h/Uio) 


-0 0940+0 0330 
U.UZ4U_o.0326 


-0 0210+°-°278 

U.UZiU_o.0274 


xVndf 


1.5 


1.5 


bins 


275 



is in a good agreement with that, extracted from the analysis of our data [|T|: A*^ = 
0.0293 ± 0.0015, i.e our data do not contradict the fi — e universality. 

In the second column the results of the joined fit of our K^s and i^'^a data are presented (this is 
practically equivalent to fixing the A+ to it's K^s value). This fit, of course, assumes the /x — e 
universality. The Aq value Aq = 0.0209 ± 0.0042 is in a good agreement with the calculations 
in the framework of the chiral perturbation theory(xPT) : Ag'^ = 0.017 ± 0.004. 
All the errors presented are from the "MINOS" procedure of the "MINUIT" program and 
are larger than the Gaussian ones. At present, we estimate an additional systematics error in 
A+, Aq to be ±0.002. The estimate is done by varying cuts, cell size during the fit of the Dalitz 
plots etc. 

In the second and in the third lines the scalar and the tensor terms are added into the fit. 
As it is seen from the second line of formula (2), the fs term is 100% anti-correlated with the 
V-A contribution (m^/2m/f)/_, i.e an independent estimate of /_ is necessary to extract fs- 

2 2 

By definition, /_ = /+(0)(Ao — A+) '"^2"^'' • A+ is, in fact, defined by the K^s data, and Aq is 
calculated by xPT- Xq^ = 0.017 ±0.004. In our fs fit we fix Aq to this, theoretical, value. The 
error(±0.004) in the theoretical prediction induces an additional error in fs equal to ±0.005. 

A possible example of theories, which give nonzero fs are the 2HDM |Ty] and the Weinberg 
3HDM model [|ll[]. In these theories, fs comes from the diagram with the charged Higgs boson 
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exchange H . The calculation of the contributions gives 
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;|-/;^(0).^.^.,/(,) (4) 

Here mn is the charged Higgs-boson mass (mass of the lightest in case of 3HDM); tg{i3) = 
V2/V1- the ratio of the vacuum expectation values for 2 Higgs doublets; a and 7 are complex 
couplings of the 3HDM Higgs boson to d-quarks and leptons. 
From our limit for fs '■ 

^ - 0.39 ± 0.2{stat) ± 0.2{theory)GeV-^; 



TTtH 



Re{al-fi)^ = 0.036 ± 0M7{stat) ± 0M7{theory) 

Hi 

Our 2HDM limit is comparable with that from LEP searches for the decay b tUj. [O: 90% 



CL. limit is < 0.4 -h 1 GeV ^(depending on collaboration). 



The results of the fit with the tensor term are presented in the third line. The tensor term 
is also correlated with Aq, dfx/dXQ = —3.5. That's why we decided to apply the same approach 
for the tensor term as for the scalar one, i.e Aq is fixed to it's theoretical value and the induced 
error in Jt, due to the theoretical error in Aq is calculated. The error equals ±0.02 for the 
single Kfj^s fit and ±0.014 for the combined one. 

The tensor coupling fx appears naturally in the leptoquark models, as a result of the Fierz 
transformation [Q. Unfortunately, we have not found complete theoretical consideration for 
this contribution. 

6 Summary and conclusions 

The decay has been studied using in-flight decays of 25 GeV K~ , detected by "ISTRA±" 
magnetic spectrometer. Due to the high statistics, adequate resolution of the detector and good 
sensitivity over all the Dalitz plot space, the measurement errors are significantly reduced as 
compared with the previous measurements. The A+ parameter of the vector formfactor 
is measured to be: 

A^ = 0.0321 ± 0.004(stat) ± 0.002(s?/st). 
and is in agreement with that obtained from our data: 

A^ = 0.0293 ± 0.0015(stat) ± 0m2{syst). 
The combined fit of both sets of data assuming the fi — e universality gives: 

A+ = 0.0296 ± 0.0014(stat) ± 0m2{syst). 
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The Aq parameter of the scalar formfactor fo{t) is measured to be: 

Ao = 0.0209 ± 0.004(stot) ± 0.002{syst). 

It is, at present, the best measurement of this parameter. It is in a good agreement with xPT 
prediction. 

The limits on the possible scalar and tensor couplings are derived: 

fs/f+{0) = 0.0039 ±0.005(stat) ±0.005(theory) ; 

fT/UiO) = -0.021 ± 0.028(stat) ±0.014(theory) 

The second(theoretical) error comes from the uncertainty in the xPT prediction for Aq. Again, 
this is the current best estimates for these parameters. 

We would like to thank V.V. Braguta, A.E. Chalov, A. A. Likhoded, A.K. Likhoded, 
R.N. Rogaliov, S.R. Slabospitsky for the discussions. 
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